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Abstract. We study the zero-temperature ground-state (GS) properties of the spin- 
i anisotropic planar pyrochlore, using the coupled cluster method (CCM) implemented 
to high orders of approximation. The system comprises a J1-J2 model on the 
checkerboard lattice, with isotropic Heisenberg interactions of strength Ji between 
all nearest-neighbour pairs of spins on the square lattice, and of strength J 2 between 
half of the next-nearest-neighbour pairs (in the checkerboard pattern). We calculate 
results for the GS energy and average local GS on-site magnetization, using various 
antiferromagnetic classical ground states as GCM model states. We also give results for 
the susceptibility of one of these states against the formation of crossed-dimer valence- 
bond crystalline (CDVBC) ordering. The complete GS phase diagram is presented 
for arbitrary values of the frustration parameter k = J2/ Ji, and when each of the 
exchange couplings can take either sign. 
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1. Introduction 

The study of frustrated quantum spin systems, both theoretically and experimentally, 
has become a held of intense activity in recent years, especially in the context of 
quantum phase transitions mEi. In particular, the experimental investigation of a wide 
variety of quasi-two-dimensional materials with fascinating properties has progressed 
hand-in-hand with the theoretical study of spin-lattice models believed to capture 
their most important behaviour. The interplay between frustration, both geometrical 
and dynamical, and quantum huctuations can produce interesting zero-temperature 
(T = 0) ground-state (GS) quantum phase transitions. These can involve phases 
ranging from ones with quasiclassical magnetic long-range order (LRO), such as Neel 
antiferromagnetic (AFM) order, to others such as valence-bond solids and spin liquids, 
which have no classical counterparts. The stable GS phases are strongly influenced by 
parameters such as the spin quantum number s of the magnetically active ions situated 
on the sites of the lattice, the dimensionality, the coordination number [z), and geometry 
of the lattice, the number of magnetic bonds and whether they are ferromagnetic (FM) 
or AFM in nature. 
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One of the key motivators that has spurred continued interest in highly frustrated 
quantum spin-lattice systems, has been the possibility of finding situations where a novel 
phase with no classical counterpart forms the stable GS phase. In particular, the search 
for quantum spin-liquid phases has received intense interest, ever since they were first 
proposed by Anderson [Ij over 40 years ago. Even though we now know that the example 
proposed by Anderson, namely, the GS phase of a spin- 1/2 Heisenberg antiferromagnet 
(HAEM) on a triangular lattice (i.e., with spins interacting via nearest-neighbour (NN) 
isotropic Heisenberg AFM exchange interactions only), is not a spin liquid, the search 
for such states continues in other systems. It is widely believed that prime candidates 
in this context are those which at the classical (s —)■ oo) level display a GS phase in 
some region of their phase diagram that has macroscopic degeneracy. In particular, 
much attention has thereby focused on frustrated quantum spin systems that are built 
from tetrahedra coupled into two-dimensional (2D) or three-dimensional (3D) lattice 
networks [51126]. 

The pyrochlore lattice comprises a 3D arrangement of corner-sharing tetrahedra, 
and it is well known that essentially all compounds that crystallize into the pyrochlore 
structure display unusual magnetic properties. In order to simplify the study of the 
3D pyrochlore, but without losing any of its magnetic frustration, it is also common to 
project the 3D vertex-sharing lattice of tetrahedra onto a 2D plane. Each tetrahedron 
has four spins at its vertices, with each of its six edges or links symbolizing a Heisenberg 
interaction. In the 2D projection (viz., the planar pyrochlore) each tetrahedron is 
mapped to a square with spins at its vertices, with the sides denoting the NN bonds with 
Heisenberg exchange bonds of coupling strength Ji, plus additional exchange bonds of 
coupling strength J 2 across its diagonals, i.e., now joining next-nearest-neighbour (NNN) 
bonds in the square-lattice geometry. Each such square is then coupled to others of the 
same kind by NN Ji bonds, resulting in the checkerboard pattern shown in figure [TJ 

The isotropic planar pyrochlore is simply the case J 2 = Ji, but it is of considerable 
interest to study also the anisotropic case when J 2 7 ^ Ji. While the 2D projection of 
the 3D pyrochlore pattern preserves the vertex-sharing structure, nevertheless one loses 
the symmetry between the six equivalent bonds on each tetrahedron in the 3D structure 
when the 2D projection is made. Thus, in the planar pyrochlore the two diagonal bonds 
of each crossed square are inequivalent to the four bonds on the sides of the square. This 
subsequent symmetry reduction is thus itself completely consistent with considering the 
case J 2 7 ^ Ji of the anisotropic Heisenberg model on the 2 D checkerboard lattice or, 
equivalently, the anisotropic planar pyrochlore. Alternate names are the crossed chain 
model and the J 1 -J 2 model on the checkerboard lattice. 

Most interest to date on the anisotropic planar pyrochlore has focussed on 
the case when both bonds are AFM in nature (i.e., Ji > 0, J 2 > 0). On 
the other hand, there has also been increasing interest at both the experimental 
and theoretical levels in quasi-2D magnetic material in which the NN coupling is 
ferromagnetic (Ji < 0), while NNN coupling is antiferromagnetic (J 2 > 0). Examples 
of such materials include Pb 2 V 0 (P 04)2 [27II3T] . (GuGl)LaNb 207 [32], SrZnV 0 (P 04)2 
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(e) 

Figure 1. The J 1 -J 2 Heisenberg model on the checkerboard lattice, showing (a) the 
ferromagnetic (FM) state, (b) the Neel state, (c) one of the two Neel* states, and (d) 
one (viz., the columnar) of the two striped states. The solid (black) lines are NN Ji 
bonds and dashed (green) lines are NNN J 2 bonds. The spins on each lattice site are 
portrayed by the (brown) arrows, (e) The classical T = 0 phase diagram. 


[SniESHSS], BaCdV 0 (P 04)2 [SSlEllES], PbZnV 0 (P 04)2 |3S] and, (CuBr)LaNb 207 [37] . 

Experimental studies of these and other materials have, in turn, stimulated theoretical 
interest in the GS and thermodynamic properties of the J 1 -J 2 model on the square 
lattice when the NN exchange bond is FM in nature (Ji < 0) and the NNN exchange 
bond is AFM, and hence frustrating, in nature (J 2 > 0) [381150] . It is of particular 
interest to note in this context that there has been a considerable degree of controversy 
between various theoretical studies on whether or not a magnetically disordered, spin- 
nematic phase emerges in the spin-| model between the phases exhibiting quasiclassical 
collinear striped AFM order and FM order. In a similar vein we have also recently 
studied ourselves the spin-| model on the honeycomb lattice with FM NN 
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bonds (of strength Ji < 0) and AFM NNN and next-next-nearest-neighbour bonds of 
equal strength (J 3 = J 2 >0) m- 

Given the high level of interest in the frustrated spin-| J 1 -J 2 Heisenberg 
ferromagnet on the square lattice (i.e., with Ji < 0, J 2 > 0) noted above, it seems 
timely to consider now the analogous model on the checkerboard lattice. In a recent 
paper [26] we studied the frustrated spin-^ HAFM on the checkerboard lattice (i.e., with 
Ji > 0, J 2 > 0), using the coupled cluster method (CCM). Our aim now is to extend 
that analysis to investigate the entire T = 0 GS phase diagram, when both exchange 
couplings Ji and J 2 can take either sign. After describing and discussing the model 
itself in section [21 we discuss the GGM methodology in section [31 Our results are then 
presented in section [H and we end with a summary and conclusions in section [5l 

2. The model 

The Hamiltonian for the anisotropic planar pyrochlore is that of a J 1 -J 2 model on a 
checkerboard lattice. It is written as 

H Ji ^ ^ Sj ■ Sj -|- J 2 ^ ^ Sj ■ Sfc , (1) 

(* 4 ) 

where the index i runs over all N sites of a square lattice; j runs over all four NN sites 
to site i that are connected to it by Ji bonds; and k runs over only two of the NNN 
sites to site i that are connected to it by J 2 bonds in the checkerboard pattern shown in 
hgure m We are interested in the thermodynamic limit [N —)■ 00 ) of an inhnite lattice. 
The sums on (q j) and ((i, k)) count each pairwise bond once and once only. Thus, the 
alternate square plaquettes of the checkerboard model have either two NNN (diagonal) 
bonds (i.e., a cross) or none, as shown in hgure [H Each site i of the lattice carries a 
particle with a spin quantum number s, dehned by a spin operator Sj = (sf, sf, sf). The 
checkerboard model is also known as the crossed chain model, which originates from its 
crossed (in our case, diagonal in hgure [1]) sets of chains, with J 2 the intrachain exchange 
coupling constant, and A the (in our case, vertical and horizontal) interchain exchange 
bonds that connect NN pairs of spins from chains running across one another. We will 
be interested here in the case s = 1/2, but it is worthwhile hrst to discuss briehy the 
classical limit (s —)■ 00 ) of the model, in order to compare with the ehects caused by 
quantum huctuations. 

For J 2 = 0 the classical GS phase is the FM state shown in hgure [I](a) for A < 0 
and the AFM Neel state shown in hgure [I](b) for Ji > 0. Both states are unique for 
specihed values of the ordering vectors. The Neel state has every row and column joined 
by Ji bonds with NN spins oriented antiparallel to one another, and hence with each 
of the diagonal crossed J 2 -chains ordered ferromagnetically with all of the spins on a 
given diagonal oriented parallel to one another, but with spins on diagonals in a given 
direction oriented antiparallel to those on the diagonals in the perpendicular direction. 
Both the Neel and the FM states will clearly remain the stable GS phases in their 
respective domains for sufficiently small values of | J 2 I. 
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The GS energy of the classical Neel state is thus = s^(—2 Ji + J 2 ). Clearly, 

for sufficiently large values of J 2 > 0 the spins on the diagonals (or crossed J 2 -chains) 
will prefer to align antiferromagnetically. It is trivial to see that there is an inhnitely 
degenerate family (IDF) of collinear AFM phases in which every diagonal J 2 -chain 
displays Neel AFM ordering, but where every diagonal (each of which is connected by 
Ji bonds to two other crossed diagonal J 2 -chains) can be arbitrarily displaced along 
its own direction or, equivalently, where all of its spins have their directions reversed. 
All of these states have the same classical energy per spin given by Efl^-p/N = —s^J 2 , 
completely independent of the exchange coupling A. 

When Ji > 0 there is thus only one classical phase transition in the model, at 
= 1, between the Neel AFM phase and the IDF of AFM phases, where k = J 2 /J 1 
is the frustration (or planar pyrochlore anisotropy) parameter. One particular member 
of this IDF is the so-called Neel* state [26], (one of which is) shown in hgure □(c). It 
has doubled AFM order, ■ ■ • ''') along each row and column of sites joined 

by Ji bonds, such that the single-site t or j, spin in the Neel state is replaced by the 
two site or H unit in the Neel* state. The Neel* state is actually doubly degenerate 
(for a given direction of the Neel vector), since a translation along a lattice (square- 
plaquette) diagonal [i.e., a translation of one lattice vector in the horizontal direction 
plus one lattice vector in the vertical direction of hgure [11(c)] transforms one of the 
pair into the other. Two other simple members of the IDF of classical GS phases for 
K > 1 (and A > 0) are the so-called columnar and row striped states [26|, the former 
of which is shown in hgure [11(d). These have, respectively, columns or rows of spins 
joined by A bonds ordered in a FM fashion, with spins on alternating columns or rows 
oriented antiparallel to one another, and hence with Neel AFM ordering along rows 
or columns, respectively, for the two cases. Like the Neel* state, the striped states are 
doubly degenerate (for a given Neel vector), since the row and columnar states transform 
into one another under interchange of rows and columns. 

The GS energy of the classical FM state is Epj^/N = s^(2A + A), and hence there 
is a second classical phase transition at «;]][ = —1 between the FM state and the IDF of 
AFM states, when A < 0. Finally, there is a third classical phase transition between 
the Neel AFM state and the FM state at A = 0, when A < 0 (or, equivalently at 
k]]) = ± 00 , with A < 0). In summary we have three classical GS phases, namely (a) the 
Neel AFM state in the regime A > 0 and — 00 < A < A, (b) the IDF of AFM states 
for A > IA I, and (c) the FM state for A < 0 and —cxo < A < —A, all as shown in 
hgure [T](e). 

Clearly the model has several interesting special limiting cases. For example, in the 
case of AFM NN couplings (A > 0), the model reduces to the isotropic Heisenberg 
antiferromagnet (HAFM) on the square lattice as k —?■ 0, and to decoupled one¬ 
dimensional (ID) crossed isotropic HAFM chains as k ^ 00 . In between, at k = 1, 
the model is simply the isotropic HAFM on the checkerboard lattice, which is itself a 
2D analogue of the 3D isotropic pyrochlore HAFM. Compared to the classical version 
(s —?■ cxd) of the model, for which the full GS phase diagram in the A^A plane is known. 
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the s = 1/2 case is basically only well established at precisely the above three points 
when Ji > 0. 

Thus, for the spin-| HAFM on the square lattice (k = 0) there is a general 
consensus that the classical Neel AFM LRO is weakened but not destroyed by quantum 
fluctuations, such that the sublattice magnetization is reduced to about 61.5% of its 
classical value of 0.5 for the s = 1/2 case [521 - I56] . The excited states are also well 
established to be gapless integer-spin magnons for the s = 1/2 model. One then expects, 
by continuity, that this (partial) Neel AFM order will persist as the frustrating NNN 
bonds (with strength J 2 > 0) are turned on. As J 2 increases (with A > 0 held fixed) 
we expect that this order will ultimately vanish at some critical value k = Kcj, where 
the Neel staggered magnetization becomes zero. In the opposite limit, k ^ 00 (i.e., 
J 2 —00 with Ji > 0), where we have decoupled ID isotropic HAFM chains, the model 
is exactly soluble for the spin-| case. The GS phase is a Luttinger spin liquid, in which 
the quantum fluctuations have completely destroyed the classical Neel LRO, such that 
the Neel staggered magnetization vanishes. The excited states are gapless, deconfined, 
spin-1 spinons. 

In between these two limits, at k = 1, most studies now concur (and see, e.g., 
[l2l[l5l[l8l|2ni[23H26]) that the GS phase of the s = 1/2 isotropic HAFM on the 
checkerboard lattice is a plaquette valence-bond crystal (PVBG) with quadrumer LRO 
on isolated spin-singlet square plaquettes, on top of which there is an excitation spectrum 
of gapped, confined, integer-spin spinons. Again, by continuity, one expects that this 
PVBG order will persist over a (finite) range of the anisotropy parameter k on either 
side of the isotropic value k, = 1. 

On the low-K side it is still not completely settled whether there is a direct (i.e., 
hrst-order in the Landau-Ginzburg scenario) phase transition between the states with 
Neel AFM order and PVBG order at k = or, if not, there exists an intermediate 
coexistence phase with both types of ordering and, hence, with two different order 
parameters. On the high-K side we expect that PVBG order will persist out to some 
critical value Kc 2 > L The situation for k > Kc 2 (with A > 0) has been, up till 
now, the most unsettled part of the phase diagram. Various scenarios have been put 
forward, as have been summarized in our own earlier paper [26]. In particular, it is 
an obvious question to ask if any of the IDF of classical (s —?■ 00 ) AFM states that 
exist for k > = 1 , survives the quantum fluctuations present in the s = 1/2 model. 

Furthermore, if the answer is affirmative, one may then enquire as to whether the 
classical degeneracy might be lifted by the well-known “order-by-disorder” mechanism 

[HZllHH]. 

In a recent previous paper [26| we studied the s = 1/2 model on the anisotropic 
checkerboard lattice, using the GGM, in the case when both NN and NNN exchange 
couplings are AFM in nature (i.e., with A > 0, A > 0). We used various AFM classical 
ground states as GGM reference states, including the Neel state of figure [11(b), the 
NeeP state of figure [11(c), and the (columnar) striped state of figure [11(d). We thereby 
presented results for the GS energy and average local on-site magnetization (i.e., the 
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magnetic order parameter) of these states, including their susceptibilities against the 
formation of both PVBC order and crossed-dimer valence bond crystal (CDVBC) order. 
The CDVBC state is one with dimer LRO on isolated spin-singlet dimers arranged in 
a pattern of crossed dimers on alternating square plaquettes (i.e., on every second row 
and every second column). Our main findings are summarized as follows. 

Firstly, we showed that the AFM quasiclassical state with Neel ordering is indeed 
the GS phase for k < Kci ~ 0.80 ± 0.01. Secondly, we showed that although quantum 
fluctuations do indeed lift the classical degeneracy of the IDF of AFM states that form 
the GS phase for the classical version (s —)> oo) of the model for k > 1, with the Neel* 
states having a lower energy than the striped states, all of these states are actually 
magnetically unstable in the sense that their magnetic order parameters are zero (or 
negative) within very small numerical errors. Thirdly, we showed that our calculations 
preferred instead a PVBC-ordered phase for < n < Kc 2 ~ 1.22 ± 0.02, and a 
CDVBC-ordered phase for all values of k > Kc 2 - Lastly, our calculations indicated that 
both transitions (i.e., from Neel to PVBC, and from PVBC to CDVBC) are probably 
direct ones. Nevertheless, we could not entirely rule out regions of coexistence in both 
cases, although we showed that if they do exist they must be very narrow, being confined 
respectively to 0.79 < k < 0.81 and 1.20 < k < 1.22. 

The main purpose of the present paper is now to extend the above analysis to the 
entire phase space, where both NN and NNN exchange bonds can now be either AFM or 
FM in nature. A particular aim will be to examine the phase boundary of the CDVBC 
state when the NN bonds are allowed to become FM (i.e, with Ji < 0). 

3. The coupled cluster method 

The CCM is a universal method of ab initio quantum many-body theory, which is 
open to systematic step-by-step improvements via various well-defined approximation 
schemes. It has been applied with considerable success to a wide spectrum of both hnite 
and extended physical systems dehned either in a spatial continuum or on a regular 
discrete lattice [591165] . These range from atoms and molecules of interest in quantum 
chemistry to the electron gas; from closed- and open-shell atomic nuclei to infinite 
nuclear matter; from various strongly interacting quantum field theories to models 
of interest in quantum optics, quantum electronics and solid-state optoelectronics; as 
well as to many condensed-matter systems, including highly-frustrated and strongly- 
correlated spin-lattice systems of the kind considered here. 

The CCM is particularly appropriate for such quantum magnets, for which many 
of the alternative methods have serious drawbacks. For example, quantum Monte Carlo 
(QMC) techniques are hindered by the well-known minus-sign problem, which is usually 
unavoidable for frustrated spin-lattice problems. Similarly, the exact diagonalization 
(ED) method is usually confined to such relatively small finite-lattice clusters that the 
more subtle orderings present in the GS phase diagram may be difficult to detect and 
calculate accurately. The CCM suffers from neither of these problems, and has been 
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applied with great success in recent years to many different spin-lattice models of topical 
interest in quantum magnetism (and see, e.g., [26lH71l^ l56 |l6^l65] - [73] and references 
cited therein). 

Since the CCM methodology has been described in detail elsewhere (and see, 
e.g., [561EIH651E711701172]), we focus only on its key elements here. We note that the 
method is size-extensive, and can hence provide results in the thermodynamic limit 
{N —)■ oo) from the outset. The first step is always to choose a suitable, normalized N- 
body model (or reference) state, |$), on top of which the quantum fluctuations present 
in the exact GS wave function IT) of the phase under study can be incorporated at the 
next stage, described below. The conditions that |<h) must satisfy are also described 
below. The exact GS ket- and bra-state wave functions, which satisfy the respective 
Schrodinger equations, 

if|T) = h;|T); {^\H = E{^\, (2) 

are chosen to have normalizations such that (T|T) = (*h|T) = (T|T) = 1 . 

A key element of the GGM is then to parametrize these exact states in terms of the 
chosen model state via the distinctive exponentiated forms, 

|T)=e^|<h); (T| = (<h|Ae-^. (3) 

In turn, the two correlation operators, S {S), are themselves formally decomposed in 
terms of a mutually commuting set of Wbody creation operators {Cf } (annihilation 
operators {Cj' = (G/)^}), as 

S = ^5,C+; S=1 + ^5,C,-, (4) 

7^0 7^0 

where we define Cq = 1 to be the identity operator, and where the set index I denotes 
a complete set of single-particle configurations for all N particles. The model state |$) 
and the operators {G/} must be selected so that |<h) is a fiducial (or cyclic) vector with 
respect to the set of mutually commuting creation operators {G/}. The model state 
thus plays the role of a generalized vacuum state, with (T|G/ = 0 = G7l*h) ,V/ 7 ^ 0. 
Furthermore the set of states {G/|<F)} spans the entire Hilbert space to which |T) 
belongs. 

For applications to spin-lattice problems it is convenient to consider each lattice 
site k in the chosen model state |<F) to be equivalent to all others. A simple way to do 
so is to rotate (passively) each spin on each site k separately in such a way that every 
spin points downward (say along the negative ^-direction) in its own local frame of spin 
axes. Such choices of local spin coordinates leave the basic SU(2) spin commutation 
relations unchanged. However, all independent-spin product model states thereby take 
the universal form |$) = | • i), and the iV-body creation operators G/ similarly 

take a universal product form, G/ = terms of the single-spin raising 

operators, = s^ + is\. If each site carries a spin with spin quantum number s, no 
site-index kj in the product form for Cf may appear more than 2 s times (i.e., for the 
present s = | case, each index may appear at most once). Glearly the set index / thus 
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simply becomes / = {/ci, ^ 2 , • • •, n = 1, 2, • • •, 2sN}. With the choice of local spin 
coordinates thus made separately for each model state, the Hamiltonian simply needs 
to be suitably rewritten in terms of the selected spin-coordinate frames. 

The CCM correlation coefficients which completely determine all GS 

properties, are now themselves calculated by minimization of the GS energy expectation 
value functional, 

H = H[Si, Si] = {^Se-^He^\^) , (5) 

with respect to each of the coefficients {Si, Si ]\fl ^ 0}. Variation of equation (j5]) with 
respect to Si from equation (j4]) yields the coupled set of non-linear equations, 

{^Cje-^He^\^) = 0 , VJ ^ 0 , (6) 

for the set of creation correlation coefficients {iS/}. Similarly, variation of equation ([5]) 
with respect to Si from equation (|1]) yields the corresponding set of linear equations, 

m~Se-^[H,Ci]e^\^), VJ ^ 0, (7) 

for the set of annihilation correlation coefficients {5/}, once equation ([6]) has hrst been 
solved for the coefficients {5/}. Since [S, G/] = 0 by construction, we note that equation 
(171) may be re-expressed in the form of a set of generalized linear eigenvalue equations, 

mS{e-^He^ - E)C+\<!>) , VJ ^ 0 , (8) 

for the coefficients {iS/}. 

The GS energy eigenvalue E, is then just the value of H at the minimum determined 
from solving equations (jb]) and ((7]), namely 

E = ($|e-^i7e^|<h) = ($|i7e^|$). (9) 

Similarly, one may hnd the GS expectation value of any other operator in terms of the 
coefficients {Si, Si}. For example, we may calculate the GS magnetic order parameter 
M, which is dehned to be the average local on-site magnetization, 

1 ^ 

M==--{4.|Se-“5^i,rV|4.), (10) 

(=1 

where in equation (1T0|) the spins are referred to their local rotated frames. 

We note that no approximations have yet been made. However, equations ((6]) 
for the coefficients {>5/} are explicitly nonlinear, and one may wonder if truncations 
of the exponential terms are needed in practice. We note that these appear, in both 
equations (j6]) and ([8]), only in the form of the similarity transformation e~^He^ of the 
Hamiltonian. This combination may be expanded in terms of the well-known sum of 
nested commutators. Another key feature of the CCM is that this otherwise inhnite 
sum will now terminate exactly with the double commutator term. This is due to the 
basic SU(2) spin commutation relations, together with the fact that all of the terms in 
equation (jl]) that comprise the decomposition for the operator S both commute with 
one another and are simple products of single-spin spin-raising operators. Thus, all 
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terms in the expansion for e~^He^ are linked, and the Goldstone linked-cluster theorem 
is exactly preserved even when the expansion of equation (jl]) for S is truncated in 
any conceivable way. In turn, this guarantees that the CCM is size-extensive at any 
such level of truncation, so that we may work in the thermodynamic limit {N —)■ oo) 
from the very beginning. Similar considerations also apply to the evaluation of the GS 
expectation value of any operator, such as that for the magnetic order parameter M in 
equation flT0|) . 

In practice, therefore, the sole approximation made is to restrict the set of indices 
{/} retained in the expansions of equation (|4]) for the CCM correlation operators {S, ^}. 
As in our earlier work [26] on this model, and in many other applications too, we 
use here the well-tested localized (lattice-animal-based subsystem) LSUBm truncation 
scheme [26 |l47l[56l[6^[67ll73] . in which at the mth level of approximation we retain only 
the multispin-flip conhgurations {/} in equation (01) that are dehned by at most m 
contiguous lattice sites. A cluster conhguration is dehned to be contiguous in this sense 
if every site in the cluster is adjacent (or connected) in the NN sense (in the selected 
geometry) to at least one other site in the cluster. 

The number of such fundamental LSUBm conhgurations, Nf = Nf{m), may be 
reduced by incorporating the space- and point-group symmetries of the Hamiltonian and 
model state being used, as well as any conservation laws that similarly pertain to both. 
Nevertheless, the number Nf{m) is a rapidly increasing function of the truncation index 
m, and hence the need soon arises to utilize massive parallelization plus supercomputing 
resources for the highest-order calculations we undertake [67101|. For the present 
checkerboard model we employ both the Neel and Neel* states shown in hgures [11(b) 
and [11(c), respectively, as our CCM model states. We also use the basic checkerboard 
geometry to dehne the LSUBm conhgurations, in the sense that we treat both the pairs 
of sites connected by Ji bonds as well as those connected by J 2 bonds as being contiguous 
sites. For both model states we are able to perform LSUBm calculations with values of 
the truncation index m < 10. 

It is worth noting that if we were to use the basic square-lattice geometry 
(i.e., with LSUBm contiguous sites dehned only by those joined by Ji bonds), the 
number Nf{m) of fundamental conhgurations would be considerably smaller than in the 
checkerboard geometry, at the same LSUBm level. In turn, this could possibly permit 
us to perform higher-order LSUBm calculations with the same level of computational 
resource. However, this possible advantage is accompanied by the severe drawback that 
the LSUBm sequences of approximations for both E/N and M now display a marked 
staggering ehect in m = 2k, depending on whether the index k is now odd or even. The 
reason for such a staggering behaviour is clearly due to the fact that the full LSUBm 
sequence using the square-lattice geometry does not properly rehect the checkerboard- 
lattice symmetries. Such odd-even staggering effects in index m have been observed 
previously in CCM LSUBm approximations for simple (dynamically unfrustrated) 
models ca. It always has the undesirable consequence of making extrapolations (for 
the full sequence) to the exact m ^ 00 limit much more difficult and much less robust. 
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Thus, finally, the only remaining step (and the only approximation made in the 
whole CCM procedure) is to extrapolate our LSUBm sequences of approximants for a 
given GS parameter to the limit m ^ oo where the method is exact. For the GS energy 
we use the usual, well-tested extrapolation scheme [26llT71[56l[65l[69lf73l[75] . 

E{m) /N = ao + aim~'^ + a 2 m~^ . (11) 

For highly frustrated magnetic systems, particularly ones which are close to a quantum 
critical point (QGP) or for which the magnetic order parameter M is close to zero, the 
extrapolation scheme which has been found to be appropriate in many earlier studies 
(see, e.g., [26llT71[56ll691l71ff73] i is one with a leading exponent 1/m^/^, 

M{m) = Co -1- . (12) 

We note finally that, without making any truncation (e.g., by using the LSUBm 
scheme utilized here), the solutions to equations ([6]) and ([8]) are formally exact, and hence 
independent, in principle, of the reference state |<h) used (at least as long, for example, 
as 1$) shares the same quantum numbers as the exact state |T)). The only caveat is 
that this assumes intrinsically that a solution (to the full, untruncated GGM equations) 
exists for a given choice of model state |<F). In practice, of course, approximations 
must be made, e.g., by using the LSUBm scheme advocated above. Then, of course, the 
solutions at a given LSUBm level of truncation for |T) may well, in principle, depend on 
the reference state |<h). As discussed more fully below in section 01 the range of values of 
the Hamiltonian parameters for which the truncated equations have a solution usually 
does, in practice, depend both on the choice of |$) and the order m of the truncation. 
So long as the scaling procedure to the m —)■ cx) limit would be exact, the dependence 
on |<F) would then be removed. In practice, since the scaling laws are empirical, the best 
we can expect is that any remaining dependence on any states |<F) for which solutions 
to the GGM equations exist in the same parameter range will be only very weak. 

4. Results 

We now present our GGM results separately for the two cases when the model state is 
chosen to be the Neel state or the Neel* state. 

4.1. Neel model state; Neel phase boundaries 

In our previous paper on the spin-| anisotropic planar pyrochlore [26] we set Ji = -|-1 
to set the energy scale for the case of AFM NN bonds, and employed the Neel state as 
our GGM model state in order to investigate the QGP at which Neel order vanishes as 
we turn on and strengthen NNN bonds with a frustrating exchange coupling of strength 
J 2 = kJi > 0. In particular, we calculated the magnetic order parameter M at various 
LSUBm levels of approximation for m < 10, for values of the frustration parameter 
K in the range 0 < k < K*{m). Here, K*{m) is the highest value for k for which, 
at a given LSUBm level of approximation, a real solution to the GGM equations (Ej) 
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Figure 2. (a) Ground-state magnetic order parameter, M, for the Neel state of the 
spin-i J 1 --J 2 checkerboard model (with Ji = — 1, J 2 < 0) as a function of I/J 2 = — 1/k. 
The COM results for various LSUBm approximations with m = 4,6,8,10 are shown, 
(b) The scaling of the LSUBm termination points k* (m) in figure [UJa) as a function 
of 1/m, is shown according to equation (fT3l) . 


could be found. Such termination points are by now well understood. They have been 
shown [171[65l[70] to be direct manifestations of the associated QCP at which the order 
of the model state melts in the physical system under study. As is typically the case, 
we found that the corresponding values of K^{m) are greater than the associated critical 
value Kci (~ 0.80 ± 0.01), and that they approach Kc^ monotonically (as a function of 
m) from above. 

We now turn our attention hrst to the case when the NNN bonds become FM in 
nature (i.e., J 2 < 0), still keeping Ji = +1. Clearly, the two bonds no longer frustrate 
one another and we fully expect the Neel order to be preserved in such a way that the 
order parameter M monotonically increases as | J 2 I increases, with M —)■ | as J 2 ^ — cxo. 
Then, as the strength Ji of the NN bonds changes sign in the same limit J 2 — 00 , 
the stable GS phase will clearly become the FM phase via a hrst-order transition. Since 
the FM states (of both the entire lattice and the decoupled crossed ID chains in this 
limit) are exact eigenstates of the Hamiltonian of equation ([1]) for all values of the spin 
quantum number s, the phase boundary between the Neel and FM phases is expected 
to be the same for the quantum case as in hgure [I](e) for the classical (s —)■ 00 ) case. 

Nevertheless, it is of interest to see how our CCM LSUBm solutions based on the 
Neel state as model state actually conform to the above expectation. Accordingly, we 
show in hgure [2](a) the LSUBm results for the magnetic order parameter M, for the 
cases m = {4,6,8,10}, with A = —1, J 2 < 0. We see once again very clearly that 
the LSUBm results based on the Neel state terminate slightly into the FM regime, 
with the corresponding lower termination points monotonically approaching the Neel- 
FM boundary at = cxo as m —)■ cx). The actual values of the inverse of the lower 
LSUBm termination points, l/fi;*(m), are +0.045, +0.0182, +0.00948 and +0.00500 
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for m = 4,6, 8, and 10, respectively. Clearly, any reasonable extrapolation scheme is 
compatible with 1/k^{oo) = 0. For example, the scheme 

l/fi:*(m) = ao + + 02?^“^ , (13) 

shown in £gure[2](b), yields the value ao = —0.0011 ± 0.0019, where the error is simply 
that associated with the £t to the four values m = {4, 6, 8,10}. 

The fact that we hnd solutions of the CCM LSUBm equations based on the Neel 
state as reference state (with a non-vanishing value of the Neel order parameter) for 
small negative values of (i.e., intruding slightly into the FM regime), is related to 
the fact that our solutions are thus intrinsically biased towards the stability of the Neel 
phase. As one proceeds to higher LSUBm orders (i.e., as the index m is increased) 
the region of unphysical intrusion decreases, and vanishes in the m —)■ cxo limit, as 
may clearly be observed in hgure [2](a). This is a completely general feature of CCM 
LSUBm solutions, which has been observed in many previous CCM studies of quantum 
magnets. In general, if the stable GS phase (say, B) in the physical intrusion region 
based on a given CCM reference state (say. A) is itself amenable to a CCM solution 
based on another (typically classical) reference state B, then we would also typically hnd 
a similar region of unphysical intrusion of the latter B-phase solution into the physical 
A-phase region. However, in the present case, since the FM wave function illustrated in 
figure [11(a) is itself an exact eigenstate of our Hamiltonian of equation ([1]), it cannot be 
used as a CCM reference state, except in the trivial sense that all (exact or approximate) 
solutions based on it simply yield vanishing values of all of the correlation coefficients 
retained in equation (jlj). 

4-2. Neel* model state; CDVBC phase boundaries 

We now turn to our corresponding LSUBm results based on the Neel* state of figure 
[T](c) as the CCM model state. In our previous paper [26] we investigated the case of 
AFM NN bonds (Ji = +1), and we now similarly consider the case of FM NN bonds 
(Ji = — 1). In both cases we investigate the effect of frustrating AFM NNN bonds 

{J2 > 0 ). 

Results for the ground-state energy per spin, E/N, for both cases (Ji = —1 and 
Ji = -|-1), are shown in hgurelHl Figure[3](a) shows the “raw” LSUBm results for values 
of the truncation index m = 4, 6, 8,10, while hgure [3](b) shows extrapolated LSUBcxo 
results using equation fITT]) . We also show the exact result for the FM phase for the 
present s = | model, Epu/N = ;j(2Ji -|- J 2 ), with Ji = —1. In hgure [3] the curves 
without symbols attached refer to our present case (Ji = —1), while the corresponding 
curves with symbols attached refer to the case Ji = -|-1. They are taken from our earlier 
work [26], and are shown for the sake of comparison. 

Figure |3](a) clearly shows that the LSUBm results for the GS energy converge 
extremely rapidly in both cases as m is increased. Furthermore, we see that for all values 
J 2 ^ 1.5 the results for both cases are remarkably similar. Nevertheless, important 
qualitative diherences emerge for small values of J 2 , in the vicinity of the corresponding 
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Figure 3. CCM results for the GS energy per spin, E/N, for the spin-i J 1 -J 2 
Heisenberg model on the checkerboard lattice, (a) LSUBm results with m = 4,6,8,10 
based on the Neel* state. Results without symbols attached refer to the case Ji = — 1, 
and are compared to the corresponding results for the case Ji = +1, shown with 
symbols attached. The GS energy per spin, E/N = (J 2 — 2)/4, for the FM phase is 
also shown, (b) The extrapolated results using equation (ITTl) : LSUBoo(l) uses the 
data sets m = {4,6, 8,10} and LSUBoo(2) uses the data sets m = {2,4, 6,8}. 


classical phase transitions at J 2 = 1 in both cases, as shown in hgure[I](e). In hgure 
n a) both sets of results are shown out to their respective (approximately determined) 
termination points, \K^{m)\. It is evident that the value of \K^{m)\ for a given value 
of m is always less in the case Ji = — 1 than in the case Ji = +1. As we showed 
previously [26], in the case Ji = +1, the values K*(m) —)■ ^^2 ~ 1-22 ± 0.02 as m —)■ cxo. 

By contrast, the respective LSUBm values of K^{m) based on the Neel* state in the 
case Ji = —1 are —0.56, k^{6) ~ —0.60, and k^{8) ~ —0.63. These termination 

points are themselves very close to the corresponding values where the GS energy curves 
cross that of the FM state, k®(4) —0.692, k^{Q) ^ —0.677, and k®(8 ) ~ —0.671. Once 

again, as the LSUBm truncation index m is increased to inhnity, the size of the region 
in which our solutions extend into the FM regime diminishes to zero, and we expect 
that fi:*(oo) = k,^{oo) = 

We note, however, that it is computationally very costly to determine the values 
K^{m) with high accuracy, since the CCM LSUBm solutions require increasingly greater 
amounts of computational resource, for a specified level of accuracy, the closer a 
termination point is approached. For this reason the values K^{m) are appreciably more 
accurate than the corresponding values K^{m). We also note that, for the same reason, 
we have been unable to track the LSUBIO solution based on the Neel* state, in the 
case Ji = —1, down to values sufficiently close to fi;*(10) for the respective energy curve 
to have crossed the FM energy curve. A simple extrapolation of the energy crossing 
points K^{m), using a scheme of the form of equation fllll) and the LSUBm values with 
m = {4, 6, 8} yields a value Kcs ~ —0.661 for the QCP. A similar extrapolation using 
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the data set m = {3, 5, 7} gives the value Kcs ~ —0.662. 

While it is certainly true that three-parameter fits to only three m-value points, 
as quoted above, are intrinsically dangerous, the energy curves themselves are very 
smooth, even near the critical point Kc^, as one can see clearly from hgure [3l Hence, 
we may have considerable confidence in the quoted values. However, to verify this 
point we may also use the extrapolated curve based on the Neel* state with the four- 
point data set m = {4,6,8,10}. Since the calculated LSUBIO curve now terminates 
before the crossing point with the FM energy curve, the resulting extrapolated Neel* 
energy curve must itself be extrapolated the short distance to the crossing point, which 
introduces an additional error. Nevertheless, since the curve itself is again very smooth, 
and indeed almost linear over the entire range of J 2 values shown in figure [3}b), such 
an extrapolation (e.g., using simple polynomial hts) is, in fact, rather robust (e.g., with 
respect to the order of the htting polynomial used). The resulting values of for hts 
using polynomials between third and sixth orders, for example, all he in the range -0.67 
to -0.68. They are thus compatible with the (more accurate) values obtained above from 
the two hts using three m-values, which require no such extrapolation to the crossing 
point. Such considerations, together with a more detailed error analysis, leads us to our 
CCM estimate Kcs = —0.66(1). 

We observe from hgure |3](b) that the energy curves rapidly approach the J 2 —)■ 00 
limit of uncoupled spin-| ID HAFM crossed chains, even for values of J 2 > 2. For 
example, the extrapolated value of the energy per spin in the regime 4 < J 2 < 5 is 
E/N —O. 442 OJ 2 , which already is in excellent agreement with the exact asymptotic 
result, E/N = (| — ln2) J 2 ^ —0.4431 J 2 from the Bethe ansatz solution [761177] . Indeed, 
if we were to include this curve in hgure [3](b) it would lie virtually on top of the 
extrapolated LSUBcx) curves, over the whole of the range 1 < J 2 < 5, and deviating 
only very slightly from them for values J 2 ^ 1- 

Although we have been able to use the Neel* state completely successfully as a CCM 
model state, we now need to investigate the stability of the quasiclassical magnetic LRO 
in the regimes in which we have used it. To that end we now show in hgure |l](a) the 
corresponding results for the GS magnetic order parameter, M, to those shown in hgure 
m a) for the GS energy per spin, E/N. It is clear that the CCM LSUBm results for M 
converge considerably more slowly as a function of the truncation parameter m than 
do those for E/N, as is to be expected by a comparison of the respective extrapolation 
schemes of equations flT^ and flTTll . Once again, just like the previous results for E/N 
in hgure [3](a), the corresponding results for M in hgure |l](a) are seen to be remarkably 
similar for the two cases Ji = — 1 and Ji = -|-1 in the region J 2 > 2. The divergence 
between the two cases for values J 2 < 2 is completely consistent with our previous 
discussion for the GS energy results. 

In our earlier work [26] for the case Ji = -|-1 we showed that the rapid rise and then 
precipitous fall in each of the LSUBm results for M, as J 2 approaches the corresponding 
termination point for that solution, is a clear marker of the associated QCP at Kc 2 ~ 1.22. 
Similarly, the analogous behaviour seen in hgure |4](a) for the LSUBm results for M for 
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Figure 4. CCM results for the GS magnetic order parameter, M, for the spin-i J 1 -J 2 
Heisenberg model on the checkerboard lattice, (a) LSUBm results with m = 4, 6,8,10 
based on the Neel* state. Results without symbols attached refer to the case Ji = — 1, 
and are compared to the corresponding results for the case Ji = +1, shown with 
symbols attached, (b) The corresponding extrapolated results using equation (fl^ and 
the data sets m = {4,6,8,10}. 


the present case Ji = —1 is evidently associated with the QCP at Kc^ ~ —0.66 that we 
have found above in the corresponding results for the GS energy. 

In £gure|l}b) we show the respective extrapolated LSUBcx) results for M for the two 
cases, as obtained from the scheme of equation (IT^ . used with the respective LSUBm 
data sets with m = {4, 6, 8,10}. We see that M is either zero (or very close to zero) or 
negative over the entire range shown in figure IH in both cases Ji = —1 and Ji = +1. 
Our results are, in particular, completely compatible with M being zero in both cases 
in the asymptotic limit J 2 ^ 00 . This is precisely the exact result for this asymptotic 
(Luttinger spin-liquid) limit of decoupled ID HAFM chains. We thus conclude that the 
Neel* state is unlikely to be the stable GS phase for any values of the parameters Ji 
and J 2 for which we have, nevertheless, successfully been able to use it as a GGM model 
state. 

Using techniques that combine renormalization group ideas, plus ID bosonization 
and current algebra methods, together with a careful analysis of the relevant terms 
near the delicate Luttinger liquid fixed point of the ID HAFM spin chain, Starykh 
et al. [20] showed that in the large-K limit the stable GS phase exhibits spontaneous 
dimerization. The GS phase in this limit has GDVBG order with twofold spontaneous 
symmetry breaking, which comprises a staggered ordering of dimers along the crossed 
J 2 chains (i.e., along the diagonals in figure H]). This finding of a GDVBG state with 
no magnetic order was then independently conhrmed [221123] in an analysis using a SE 
technique based on the flow equation method. 

It is thus of interest to investigate within our present GGM analysis whether the 
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Figure 5. (a) CCM results for the scaled inverse crossed-dimer susceptibility, n/xd, 
as a function of the frustration parameter, k = J 2 IJi, using the Neel* state as the 
CCM model state, for the spin-i J 1 -J 2 Heisenberg model on the checkerboard lattice 
with Ji = — 1 (left curves), comparing with those for Ji = +\ (right curves) (with 
J 2 > 0 for both cases). We also show the extrapolated LSUBoo results based on 
equation (1161) with the data sets m = {4,6,8,10} as input, (b) The perturbations 
(fields) F = SOd for the dimer susceptibility Xd- Thick (red) dashed and thin (blue) 
dashed lines correspond respectively to strengthened and weakened NNN exchange 
couplings, where Od = J2{{i k)) ■ Sfc, and the sum runs over the NNN diagonal 

bonds of the checkerboard lattice, with Qik = +1 and —1 for thick (red) dashed and 
thin (blue) dashed lines respectively. The original solid (black) Ji bonds are unaltered 
in strength. 


whole (or part) of the regime that has been accessible to us using the Neel* state as 
model state, but for which we have shown has no Neel* magnetic LRO, might in fact 
have CDVBC order instead. To that end we consider the response of the system when a 
held operator, F = 60d, is added as a small perturbation to the original Hamiltonian of 
equation ([T]), with 6 an inhnitesimally small c-number, and the operator Od, illustrated 
in £gure[5](b), promotes the formation of CDVBC order. 

Thus, we now calculate the perturbed GS energy per spin, e{6) = E{6)/N, for the 
perturbed Hamiltonian, H{6) = H+F, at various LSUBm levels of approximation, using 
the Neel* state as our model state. We then calculate the corresponding susceptibility 
coefficient. 


d^e{d) 

86 ^ 


1 

( 5=0 


(14) 


which is a measure of the susceptibility of the system against the formation of CDVBC 
order. Clearly, Xd can be used to hnd points or regions in phase space where the 
phase corresponding to the particular CCM model state used (viz., here the Neel* 
state) becomes unstable with respect to a CDVBC-ordered state, namely when the 
extrapolated inverse susceptibility, 1/xd, goes to zero. 
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In order to extrapolate our LSUBm results to the LSUBoo limit, it has been 
found [78l|79] that the simple scheme, 

Xd^im) = xo + + X2m~'^, (15) 

gives excellent hts, except in regions where Xd^ becomes very small or zero. In the 
present case we are, of course, interested in precisely such regimes, and it is then 
preferable [72l[^[79] to use an unbiased extrapolation scheme of the form, 

Xd^im) =yo + yim-’', (16) 

where the leading exponent n is itself a htting parameter, together with the coefficients 
yo and yi. 

We show our results for Xd^ in figure [5](a). Since the GS energy scales linearly with 
J 2 in the large J 2 limit, as seen from £gure[3l it is more appropriate to show our CCM 
results in £gure|5](a) for the scaled inverse dimer susceptibility, n/xd, as a function of n. 
Figure [5](a) shows results for both the cases Ji = — 1 and Ji = +1. It is evident from 
LSUBcx) extrapolations using equation flTBl) that our results in both cases are completely 
consistent with 1/xd being zero for all values of k shown. In other words, our results 
strongly indicate that everywhere we have been able to use the Neel* state as a CCM 
model state it is actually unstable against the formation of CDVBC order, which forms 
the ordering of the true stable GS phase in these regimes. Hence, we conclude that 
at Kc 2 = 1-22 ib 0.02 there is a QCP between states with PVBC and CDVBC forms of 
order, and at Kc^ = —0.66 ± 0.01 there is a QCP between states with CDVBC and FM 
forms of order. 

5. Summary and conclusions 

In this and our previous paper [26] we have used the CCM, implemented to high orders, 
to study the spin-^ J 1 -J 2 Heisenberg model on the checkerboard lattice. Whereas in 
our earlier work, we studied the model only in the case where both NN and NNN bonds 
were AFM in nature (Ji > 0, J 2 > 0), we have now investigated the model over the 
entire J 1 -J 2 phase space. Previously we showed that the classical phase transition at 
= 1 (Ji > 0) between the Neel AFM phase and the IDF of AFM phases, is split in 
the s = I model into two transitions at = 0.80 ± 0.01 and Kc 2 = 1-22 ± 0.02. For 
K < Kci the Neel order persists, while for < k < Kc 2 the GS phase has PVBC order. 
Finally, we showed that for k > Kc 2 (with A > 0) the stable GS phase has CDVBC 
order. The transitions at both QCPs and k,c 2 are likely to be direct ones, although 
we cannot exclude very narrow coexistence regions conhned to 0.79 < k < 0.81 and 
1-20 < K < 1.22, respectively. 

We have now conhrmed that the CDVBC phase persists when the A bonds become 
FM in nature (A = —1) to a QCP at Keg = —0.66 ± 0.01, at which point the CDVBC 
phase gives way to the FM phase. This QCP may be compared with the corresponding 
transition in the classical (s —?■ 00 ) model at = —1 between the IDF of AFM states 
and the FM state. It is interesting to note that quantum fluctuations act to destabilize 
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Figure 6. The ground-state (T = 0) phase diagram of the spin-l J 1 -J 2 Heisenberg 
model on a checkerboard lattice (with n = J 2 /J 1 = tan0), showing the Neel 
antiferromagnetic phase, the plaquette valence-bond crystalline (PVBC) phase, the 
crossed-dimer valence-bond crystalline (CDVBC) phase, and the ferromagnetic (FM) 
phase. All of the transitions at Kci ~ 0.80(2), Kc2 ~ 1-22(2), and Kc^ ~ —0.66(1), each 
with J 2 > 0, appear to be direct ones. The direct first-order transition between the 
FM and Neel phases is exactly at 9^ = |7r. 


the FM ordering at a weaker level of frustration than for the classical version of the 
model. Precisely the same effect has now also been seen in a variety of other comparable 
models, such as the spin-| J 1 -J 2 model on the square lattice [H], the spin-^ J 1 -J 2 -J 3 
model (with J 3 = J 2 ) on the honeycomb lattice [80], and the spin-i J 1 -J 2 model on 
a cross-striped square lattice Finally, in the unfrustrated region where J 2 < 0, 

the QCP between the FM and Neel phases has been shown to occur at Kd = 00 , at 
exactly the same place as in the classical model (^^ 1 )) Mly expected. Our results are 
summarized in the complete GS (T = 0) phase diagram shown in figure | 6 l 

We note that there have been suggestions in the literature [39lll0l|44l|45lll9] that 
the competition between FM interactions between NN pairs of spins (Ji > 0) and 
AFM interactions between other pairs of spins in frustrated spin-i systems on the 
square lattice might result in gapless spin-liquid states with multipolar order (e.g., spin- 
nematic states) near to the FM boundary. Similar states have also been hypothesised in 
frustrated multiple cyclic spin-exchange models on the triangular lattice with FM NN 
pairwise interactions |19], either in a nonzero external magnetic held (with octupolar 
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ordering occurring) or in zero field (with quadratic or nematic ordering occurring in a 
state bordering the FM phase). Nevertheless, such states with multipolar-ordering in 
the zero-held case are deemed to be quite fragile. Indeed, a recent careful and accurate 
analysis of the spin-i FM version of the J 1 -J 2 Heisenberg model (i.e., with Ji < 0) 
on the square lattice m, which used both high-order CCM and ED techniques, found 
that if such states did exist (and no evidence at all was found for them in this study), 
they could exist only over a very small range below J 2 ~ 0.4| Ji|, where the transition 
at which FM ordering disappears was accurately determined to be = 0.394(1) |Ji|. 
A recent Schwinger boson study of the same model [81] also found no evidence for any 
such states. 

Similarly, in our present study of the checkerboard model we have found no evidence 
at all to suggest that the CDVBC phase does not extend all the way down to the FM 
phase. Again, if any such intermediate phase exists at all, it would need to be conhned 
to a very small range below J 2 ~ 0.7| Ji|, on all the evidence presented here, where the 
transition at which FM ordering disappears is = 0.66(1) |Ji|. On the other hand, the 
detection of phases with novel quantum ordering, such as spin-nematic states, is subtle, 
and the present checkerboard model may merit further investigation in this respect in 
the very narrow region just above the border of the FM phase. 

Finally, we note that there has been interest shown in frustrated ferromagnets 
with respect to the formation of multimagnon bound states under the influence of 
high magnetic fields (and see, e.g., [3^[821 - I84] 1. It might, therefore, be of interest 
to investigate the present checkerboard model further, by coupling it to an external 
magnetic held. 
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